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ABSTRACT 

Context. For direct imaging of exoplanets, a stellar coronagraph helps to remove the image of an observed bright star by attenuating 
the diffraction effects caused by the telescope aperture of diameter D. The Dual Zone Phase Mask (DZPM) coronagraph constitutes a 
promising concept since it theoretically offers a small inner working angle (IWA ~ A /D where A denotes the central wavelength of 
the spectral range A/1), good achromaticity and high starlight rejection, typically reaching a 10 6 contrast at 5 Aq/D from the star over a 
spectral bandwidth AA/A of 25% (similar to H-band). This last value proves to be encouraging for broadband imaging of young and 
warm Jupiter-like planets. 

Aims. Contrast levels higher than 10 6 are however required for the observation of older and/or less massive companions over a finite 
spectral bandwidth. An achromatization improvement of the DZPM coronagraph is therefore mandatory to reach such performance. 
Methods. In its design, the DZPM coronagraph uses a grey (or achromatic) apodization. We propose to replace it by a colored 
apodization to increase the performance of this coronagraphic system over a large spectral range. This innovative concept, called 
Colored Apodizer Phase Mask (CAPM) coronagraph, is defined with some design parameters optimized to reach the best contrast in 
the exoplanet search area. Once this done, we study the performance of the CAPM coronagraph in the presence of different errors to 
evaluate the sensitivity of our concept. 

Results. A 2.5 mag contrast gain is estimated from the performance provided by the CAPM coronagraph with respect to that of the 
DZPM coronagraph. A 2.2 ■ 10~ 8 intensity level at 5 A /D separation is then theoretically achieved with the CAPM coronagraph in 
the presence of a clear circular aperture and a 25% bandwidth. In addition, our studies show that our concept is less sensitive to low 
than high-order aberrations for a given value of rms wavefront errors. 

Key words. Instrumentation: high angular resolution - Techniques: high angular resolution - Telescopes - Methods: numerical 



1. Introduction based telescopes ([Macintosh et all 120081: ISauvage et al.l l2010t 



iGuvon et alJl2010bLlHinklev et al.ll201 ll) will compensate the ef- 



m p ' 1.1. Astronomical context fects of the atmospheric turbulence to provide high angular res- 

J> • olution images at very high Strehl ratio. A stellar coronagraph 

.,H , Recently, the direct imaging of several exoplanets have dGuvon et alJ I2006) will help to mostly suppress the diffrac- 

^ ; been made possibfe (|Kalas etal. | |2008fc i Marois etalj |2QQ8 tion effects due t0 the tele scope apertu re. In addition, post- 

b . l Lagrange et al. || 2010|) owing to their youth (tens of Myr), their proce ssing techniques (e . g. IMarois iet alJB oOO: Scarks & Ford 

. * , high mass (a few Jupiter masses) or their wide apparent dis- 2QQ2 . Marois et al 2006; Lafren iere et al. 2007; Mugnier et al. 

tance from their host star (larger than 1 arcsecond). For imaging 2008 , l V i gan etalJ 12001) will allow us to remove the speckles 

older, less massive or less separated planetary companions while caused by the XAO residual aberrations and the slow variations 

enabling their spectroscopic observations, achromatic high con- of the physical conditions of a telescope. 

trast imaging techniques are required. However, large broadband The current Hmits of all mesg methods need to be pushed 
imaging of these exoplanets is extraordinarily challenging in the further to expect imaging me ex0 planets mentioned above, 
visible and near-infrared bands because of the large flux ratio In particulari increas ing the performance of the current coro- 
(10 7 to 10 10 ) and the small angular separation (less than 1 arc- nagraphs is man datory for reaching an excellent polychro- 
second) between them and their host star ( |Traub & OppenheimeiJ matic removal of the star diffraction p attern and image the 
| /U1H ). faintest substellar mass companions in white light. The re- 
Such detection requires the combination of different methods, sulting> innovative coronagraphs will find their place in the 
such as extreme adaptive optics (XAO), stellar coronagraphy exoplanet imagers planned to follow the upC oming Palomar 
and post processing techniques, to enhance the exoplanet im- pl64Q (Hinklev et al 2 011), VLT-SPHERE (Beuzit et al. 2008), 
age (Oppenheimer & Hinklev 2009). XAO systems on ground- 
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Gemini Planet Imager ( Macintosh et alJ 120081) and Subaru 
HiCIAO (iHodann e t al. 2008): e.g. EPICS dKasneret alJl2008h 
for the European Extremely Large Telescope (E-ELT) on 
the ground or possible l-2m off-axis telescopes in spa ce 
dGalicher et al.ll2010tlTrauger et al.ll20lol; iGuvon et ai]|2010d) . 

1 .2. Broadband capabilities of coronagraphs 

Many coronagraphic designs have been proposed these past few 
years; most of them are reported in lGuvon et al.l d2006l) . We are 
here particularly interested in broadband behavior of corona- 
graphs and although a complete review of the broadband per- 
formance of all the existing coronagraphs is not within the scope 
of this paper, we provide here some numbers from the literature 
allowing to situate the presented work in the context of the state 
of the art. Aq and D denote the central wavelength of the spectral 
range AA and the telescope aperture diamet er respectively. 
The Achromatic Interfero Coronagraph (Gav & Rabbia 1996) 
was the first attempt to make a concept able to image faint sub- 
stellar mass companions close to the star in white light. This 
device can theoretically offer an intensity level of 1.6 ■ 10~ 5 
at a dist ance of Aq/D from t he star over a bandwidth (AA/Aq) 
of 18% (iBaudozet all l2Q00h. The Apodized Pupi l Lyot coro- 
nagraph (Aim e et al.ll2002t ISoummer et al.ll2003al) . a more re- 
cent concept, is one of the most popular designs, adopted for 
the forthcoming generation of exoplanet imagers (e.g. VLT- 
SPHERE, Gemini Planet Imager, Palomar P1640). It can the- 
oretically reach an intensity level lower than 10~ 7 at 5 Aq/D over 
a bandwidth of 20% , in the presence of a telescope aperture with 
centra l obstruction ( Aim el2005t ISoummerl2 005 ; Sou mmer et al.l 
2011b. Band-hmite d coronagraphs (Kuchner & Traubl 12002: 
Kuch ner et al.l2 005) also represent encouraging solutions to pro- 
duce very high contrast performance over large spectral bands. 
Recent studies predict for instance that these concepts using hy- 
brid metal-dielectric masks and combined with a deformable 
mirror, set to generate a dark half field in the image plane, 
can provide an intensity level of 4.9 • 10~ 1 at 3.5 Aq/D withi n 



this dark region, over a bandwidth of 20% ( Moo dy et al. 
The Phase Induced Amplitude Apodization (PIAA, 
2003) and its avatars (Guvon et al. 2005; Martinache et a 



20081) 



Guyo 
.1200* 



Pluzh nik et al.ll2.006t IGuvon et al.ll2010al) constitute some of the 
most promising solutions for Earth-like planet detection. Indeed, 
an intensity level of 10~ 1() at 1.5 Aq/D over a 21% bandwidth is 
predicted with an hyb rid coronagraph comp osed of PIAA with a 
classical apodization dPluzhnik et al.ll2006l) . 
Stellar coronagraphs using a phase mask in the focal plane 
represent some encouraging solutions to achieve the re- 
quested contrast levels, either i n the form of secto rized masks 
(Four Quadrant Ph ase Mask j Roiian et al.l l2000ll and Eight 
Octant Phase Mask | Murakami et al. 2008 1) or ci rcular masks 
dRoddier & Roddierl 1 1997b ISoummer et alj|2003bl) . These sys- 
tems are characterized by their very small inner working angles 
(IWAs) and good starlight rejection capabilities. 
Achromatization of these phase mask coronagraphs has been in- 
tensively studied in the past few years to improve the capability 
of these concepts to work in broad spectral bands. For sectorized 
phase masks, the use of achromatic half-wave plates leads to in- 
tensity leve ls of 5 • 10~ 5 at 2. 5 Ap/D from the star over a 20% 
bandwidth dMawet et al.ll2006h . 

Derived from the FQ PM concept, the vector vortex coronagraph 
dMawet et al.l 120051 120 1 Oh also provides achromatic solutions. 
Using Liquid Crystal Polymers (LCP) for the manufacturing of 
such phase mask, intensity levels lower than 1.4 • 1 0~ 7 at 3 Aq/D 
for a 15% bandwidth can theoretically be reached ( Mawet et alJ 



l2005l) . and fur ther improvements are expected using a 3-layer 
LCP structure dMawet et al.ll2010h . 

Achromatization of circular phase mask coronagraphs by the 
use of dual zone ph ase structures has also been studied 
dSoummer et alJ 1200 3b). leading to intensity levels of 10~ 6 at 
3A/D for 20% bandwidth. We report on further improvements 
of this concept in the following. 



1.3. Circular phase mask coronagraphs 

Focus is made here on stellar coronagraphs with circu- 
lar axi-symmetrical p hase masks dRoddier & Roddierl 1 1997; 
ISoummer et alJl2003bl) . This is a promising concept, well suited 
for any arbitrary aperture (clear aperture, centrally obstructed 
circular pupil, diluted pupil...) and free from blin d axes in the 
coronagraphic images. Rodd ier~& Roddierl d 19971) proposed the 
first design, referred here as Roddier & Roddier Phase Mask 
(RRPM), which was later improved by adding an entrance pupil 
apodizer (Guvo n~& Roddierll2000t [Soummer et ID|2003a). The 
corresponding system, called apodized pupil Roddier & Roddier 
phase mask (ARPM) coronagraph, can theoretically achieve per- 
fect starlight extinction at a single wavelength. Unfortunately, its 
performance decreases strongly in broadband observation. This 
loss of coronagraphic performance comes from chromatism ef- 
fects related to t he inhe rent properties of the RRPM. 
Sou mmer et alJ d200 3b) proposed an improvement of the ARPM 
coronagraph to overcome the problem of chromaticity, replacing 
the simple phase mask of the Roddier coronagraph with a dou- 
ble phase mask, named Dual Zone Phase Mask (DZPM). With 
this innovative coronagraph which also includes pupil apodiza- 
tion, they showed that an intensity level of ~ 10~ 6 can theo- 
retically be reached at a 3> Aq/D angular distance from the star. 
These promising results were achieved both for a clear circular 
aperture and for an entrance pupil with 14% central obscuration 
ratio, in both cases considering a 20% bandwidth. Although en- 
couraging, the performance of this coronagraphic system needs 
to be enlarged so as to realize spectro-imaging observations of 
the faintest close-in companions with this concept. The DZPM 
coronagraph has been originally designed considering a grey (or 
achromatic) entrance pupil apodization. This component then 
applies the same transmission function to the incident beam over 
the whole spectral bandwidth of study. The apodization shape is 
possibly not optimal at all the wavelengths of study. 



1.4. Objectives of the paper 

The aim of this paper is twofold. First of all, we investigate the 
possibility of using a colored (or chromatic) instead of a grey 
apodization to improve the contrast gain delivered by the DZPM 
coronagraph over the whole spectral bandwidth. This innovative 
concept will be named Colored Apodizer Phase Mask (CAPM) 
coronagraph. Secondly, once the apodization parameters opti- 
mized, we propose to analyze the sensitivity of the CAPM coro- 
nagraph by estimating the impact of different effects which may 
alter its performance in a realistic case. 

In Section [2] we briefly review the principle of the DZPM coro- 
nagraph and state the possibility of upgrading this concept with 
the use of a colored apodization. The methodology, applied to 
optimize numerically the parameters of the CAPM coronagraph, 
is thereafter described in Section [3] We present the theoretical 
performance achieved by this coronagraphic device for different 
spectral bandwidths and a clear circular aperture in Section|4] In 
addition, manufacturing aspects are addressed in Section with 
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the development of a physical design for the colored apodiza- 
tion and the performance study for the CAPM coronagraph in 
the presence of this component. The sensitivity of the CAPM 
coronagraph to different errors is thereafter analyzed in Section 
i 



2. The CAPM coronagraph 

The principle of the DZPM coronagraph is recalled in the fol- 
lowing before introducing our proposal, the CAPM coronagraph 
and its colored apodization. We then expose the formalism of 
this latter. 



2.1. Principle of the DZPM coronagraph 

An illustration of the DZPM coronagraphic layout is given in 
Figure [T] detailing the typical four planes of Lyot style corona- 
graphs. The DZPM coronagraph combines an apodization in the 
entrance pupil plane A, a DZPM in the following focal plane 
B and a Lyot stop in the re-imaged pupil plane C. The coro- 
nagraphic image forms in the final image plane D in which a 
detector is located. The optical design is such that the complex 
amplitudes of the electric field in two successive planes are re- 
lated one to each other by a Fourier transform operation. 
The DZPM is designed as a circular phase disk of diameter d\ 
surrounded by an annular phase ring of diameter d 2 , see Figure 
[2] Each diameter is equal to a very precise fraction of the size of 
the Airy disk. The inner and outer parts of the mask introduce 
some phase shifts tp\ and (f2 respectively, to the incoming wave- 
front in the focal plane B. As the wavelength increases, <p\ and 
if 2 decrease and the stellar diffraction pattern expands. This com- 
bination of effects leads to a decrease of the relative flux passing 
through the inner zone, while the flux passing through the annu- 
lar zone increases. At a given wavelength, the complex ampli- 
tude of the field in the relayed pupil plane C corresponds to the 
result of the interference between the direct wave with the waves 
delayed by the two phase shift zones of the DZPM. This sum of 
wavefronts varies with the wavelength since <p\, <p2 and the ar- 
eas of the diffraction pattern covered by the mask also evolve 
with the wavelength. With an adequate choice of the DZPM pa- 
rameters, the recombination of these three wavefronts in plane C 
generates destructive interferences inside the geometric pupil for 
a wide range of wavelengths. The light is rejected outside of this 
relayed pupil and blocked by a Lyot stop, leading to a broadband 
suppression of the star image at the detector in plane D. 
In addition, the presence of an entrance pupil apodization in 
plane A allows us to match well the direct wave with the waves 
diffracted by the DZPM. It results in improved destructive inter- 
ference over the relayed pupil. A lower amount of the residual 
starlight is then found inside this pupil and therefore, a better 
starlight extinc tion is observed in the fi nal image plane D. As 
was shown by Soumm er et al.l (l2003bl) . further improvements 
can be achieved by adding a slight defocus of the coronagraph 
mask, which is equivalent to allocate a complex term in our 
apodization. 



2.2. Replacement of the grey by a colored apodization 

A complex apodization in the coronagraphic scheme allows the 
direct wave and the waves diffracted by the DZPM to match well 
within the geometric pupil image and inte rfere destructively over 
it. In their design, Soum mer et al.l d2003bl) considered a grey ver- 
sion of the complex apodization. Since the waves diffracted by 
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Fig. 1. Scheme of the coronagraphic layout. 
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Fig. 2. Schematic representation of the DZPM. The dimensions 
are not at scale in this drawing. 



the mask evolve with the wavelength, they do not match well the 
direct wave at all the wavelengths. Ideally, a complex apodiza- 
tion has to be found for each wavelength to obtain an achromatic 
nulling recombination of the direct with the diffracted waves. 
The design improvement proposed here from the DZPM coro- 
nagraph is similar to the approach developed by lAimd ((2005) 
for the Prolate apodized Lyot coronagraph. The idea consists 
in replacing the grey by a colored apodization in our corona- 
graphic system. This colored apodization constitutes an ampli- 
tude apodization and hence, intensity apodization variable with 
the wavelength. This change in apodization leads to the design 
of CAPM coronagraph which combines a colored apodization, 
a DZPM in the following focal plane and a Lyot stop in the re- 
layed pupil plane. With this innovative coronagraphic system, 
we expect to reach contrast levels larger than those achieved by 
the DZPM coronagraph over a finite spectral bandwidth. 

2.3. Formalism 

In the following, the formalism common to DZPM and CAPM 
coronagraphs is given. For the sake of clarity, we omit the posi- 
tion vector r, its modulus r and the wavelength A in the equations 
below. T symbolizes the Fourier transform operator in which we 
include the Fourier Optics scaling factor with / the tele- 

scope focal length dGoodmanll996l) . 

The complex amplitude of the field *Fa at the aperture is given 
by: 

x Pa=P<S>, (1) 

in which P defines the telescope aperture shape whereas €> de- 
notes the complex apodization. 

The DZPM is located in the following focal plane B. Its ampli- 
tude transmission function t can be written as: 



t = 1 - (e m - e^' ) Mi - (1 - e m ) M 2 , 



(2) 
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in which Mi and M2 define the top-hat functions of the inner and 
outer parts of the DZPM; they are equal to 1 for \r\ < d\/2 and 
|r| < c/2/2 respectively and otherwise. Let us remind that tpi 
and ip2 represent the phase shifts introduced by the mask in its 
inner and outer parts respectively, at a given wavelength. 
The complex amplitude of the field after the mask is then 
given by: 

V B = [l - (e^- - )M X - (1 - e'n M 2 ] TTSjd . (3) 

The DZPM and CAPM coronagraphs belong to the class of the 
Lyot-style coronagraphs which consist in a succession of binary 
filters (pupil, finite-size focal plane mask, Lyot stop). We then 
adopt the semi-analytical approach suggested by Soummer et al.l 
(2007) in the case of a Lyot-style coronagraph for a better under- 
standing of the interference effects occurring in the pupil plane 
C. The complex amplitude after the Lyot stop *Fc is expressed as 
follows: 

- (e^ - «** ) T [TWaWi] L (4) 
-{\-e^-)T[TVY A W2\L, 

in which L denotes the index function of the Lyot stop. *Fc is 
the sum of three terms which represent the direct wave ^ and 
the sum of the waves diffracted by the inner and outer parts of 
the DZPM. As mentioned above, perfect starlight cancellation 
is achieved if the sum of the waves diffracted by the DZPM 
matches perfectly the direct wave inside the geometric pupil. 
The coronagraphic amplitude in the image plane D is finally 
obtained with: 

= fVfc\ • (5) 

The corresponding intensity Id of the polychromatic corona- 
graphic image is given by: 

Id = -t7 \ W D \ 2 dA. (6) 

A ^ Ja„-AA/2 

We recall that AA defines the spectral range of study centered at 
the wavelength Ao. 

3. Methodology 

In this section, we describe the methodology used to realize our 
numerical simulations of the CAPM coronagraph. We also ex- 
pose the optimization criteria to determine the parameters of this 
latter. 



of the central obscuration for a centrally obstructed circular aper- 
ture. For a given central obscuration ratio 77 (null for a clear aper- 
ture), the amplitude apodization can be written as: 

<D fl (r, A) = 1 + co u [r 2 - + cj 2 .a (r A - , (8) 

where co\x and o>2,a denote the colored parameters related to the 
amplitude apodization <!>„. In addition, the phase apodization <t„ 
is expressed as follows: 

<D vv (r, A) = exp(2i7i/3 A r 2 A Q /A) , (9) 

where fix is a chromatic parameter associated with the phase 
apodization. The chromatic longitudinal defocus Az is related 
to fix by the following expression: 

Az(A) = 2f3 A f 2 A, (10) 

where / denotes the optical focal length of the system. 
Manufacturing aspects on the chromatic defocus are discussed 
below in Section [5] As mentioned above, oj\j, u>2,a and fix are 
chromatic parameters since we deal here with a colored apodiza- 
tion. They will be optimized at each wavelength of the spectral 
bandwidth with the criteria given below. Let us note that the pre- 
vious expressions of apodization are still valid for the DZPM 
coronagraph, with a>ij, u>2,a and fix remaining constant over the 
spectral bandwidth. 

3.2. Numerical simulations using ultrafast Fourier transforms 

We detail here the computation of the amplitudes *Pc an d ^d, 
respectively based on a semi-analytical approach and a novel 
method providing fast computation of high resolution Fourier 
transform. 

3.2.1 . Computation of the amplitude in plane C 

The complex amplitude in pupil plane C is computed fol- 
lowing Eq.©. This semi-analytical expression is advantageous 
in terms of computation time since Fourier transform operations 
are only required within the mask area of width ~ Ao/D in the 
focal plane B and within the Lyot stop area in plane C. No zero- 
padding is required in this case so we can use arrays of width 
eq ual to the mask d iamet er. We apply the algorithm, proposed 
by ISoummer et al.l (I2007I) . based on the matrix Fourier trans- 
form (MFT) to compute these truncated Fourier transform and 
evaluate SCc- 



3.1. Parameters of the CAPM coronagraph 

Several parameters are involved in the CAPM coronagraph. d\, 
d 2 and the phase steps, expressed in optical path differences 
OPDi and OPD2, correspond to the characteristics related to the 
DZPM. The other parameters refer to the complex apodization 
<t> which can be decomposed into the product of an amplitude 
and a phase apodizations <t> fl and <t>„. respectively: 

<D = <D fl ■ O w . (7) 

Following the indications of Soumme r et alj d2003bl) . the trans- 
mission function chosen for <t>„ is a radial symmetric fourth or- 
der polynomial while <D„ is in the form of a defocus of the mask. 
This amplitude transmission function finds its maximum value 
at the pupil center in the case of a clear aperture and at the edge 



3.2.2. Computation of the amplitude in plane D 

The complex amplitude To in image plane D is obtained with 
a Fourier transform of the complex amplitude Tc, see Eq.©. 
We want to evaluate this complex amplitude over a large field of 
view (FoV), typically from 50 to 200 Aq/D width with good sam- 
pling. This can be seen as a truncated Fourier transform within 
the area of interest (here, FoV) like in the case of and the 
MFT can be used again. However, the computation proves to be 
slow since an array much larger than that used in plane B is in- 
volved. While this increase in computation time is not dramatic 
in the monochromatic case, it becomes important when we con- 
sider a broadband image for which several monochromatic im- 
ages are computed and summed. 

For fast computation of from Wc, we apply a novel method 
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based on a succession of Fast Fourier Transforms (FFTs) in- 
stead of the MFT. By this method, already experimented in 
the context of James Webb Space Telescope PSF calculation 
dZamkotsian & Doh len 2004), the diameter in pixels of the ge- 
ometric pupil in plane C is assumed to be equal to Nco- is 
represented in an array Co with dimensions Nco x Nco- In the 
following, we describe the successive steps used for the compu- 
tation of and illustrated in Figure[3j 

1) FFT of the array Co- This leads to an array D\ which contains 
*Po over a certain FoV of width Nd\ A/D where Nd\ = Nco- 

2) Selection of a central sub-array of D\ to restrict the FoV to the 
area of interest of width Nm A/D. The dimensions of this new 
array D2 are Nm X Nm with Nm ^ Nd\- It contains Td in a 
smaller FoV than in D\ . 

3) Inverse FFT of the array D^- It results in an array C2 of same 
size as D2 containing a miniature version of the pupil. 

4) Use of the zero padding technique with the insertion of C2 in 
a larger array C3 of dimensions Nc3 X Nci ■ 

5) FFT of the array C3. The result of this operation is an array 
Z>3, of the same dimensions as C3, containing the final complex 
amplitude *¥d in plane D. 

Thanks to these successive operations, a fast computation of 
is achieved. To check the accuracy of our method, we compare 
the monochromatic intensity |*To| 2 obtained with our ultrafast 
Fourier transform method and with a simple FFT (steps 1 to 3 are 
skipped here). We work with Nco = 390 and Nqj = 10 Nco so 
the element resolution A/D is sampled each 10 pixels in the im- 
age plane. This leads to an absolute difference less than 2 • 10~ 15 . 
This numerical noise remains largely below all the contrast val- 
ues that we are considering here and therefore, it does not have 
an influence on the analysis of our results. This test allows us to 
confirm the accuracy of our ultrafast Fourier transform proposed 
here. In the following, we apply this algorithm to estimate Wd- 
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Fig. 3. Schematic illustration of the method used for the compu- 
tation of 



3.3. Optimization criteria 

In their study. ISoummer et alJ d200 3b) optimized the parameters 
of the DZPM coronagraph to achieve the best on-axis extinction 
of the stellar diffraction pattern in the image plane. We decide 
here to optimize the CAPM coronagraph parameters in order to 
reach the best off-axis cancellation of the star image. This will 
yield minimal starlight pollution in the search area in which ex- 



oplanets are expected to be found. The search area is defined by 
an annulus of width between 2 Ao/D and 10 Ay/D from the main 
optical axis in the image plane. 

We average the intensity of the monochromatic coronagraphic 
image over this search area. This calculus is performed at five 
different wavelengths uniformly distributed within the spectral 
bandwidth. Once this done, we compute a merit function M to 
estimate the quadratic sum of the mean intensities reached at 
each wavelength A: 

K Po 

M^^l^l 2 , (11) 

<<i Pi 

where pj and po are the inner and outer radii of the search area 
and A\ and A„ represent the first and n-th wavelengths of study 
within the spectral bandwidth. For our optimization, 5 wave- 
lengths from A\ to As have been considered within the spectral 
bandwidth AA, centered on A3 - Aq, and equally spaced from 
one to another by AA/5, see Figure [4] Some test optimizations 
have also been run with an eleven wavelengths gridding for 25% 
and 50% bands, showing no improvement of the coronagraph 
performance with respect to the five wavelength sampling. 
A numerical least-squares method is used to obtain the CAPM 
coronagraph parameters (d\, d.2, <pi, <fi2, <^>\,a, ol>2,a, Pa) that min- 
imize our merit function and therefore, achieve the lowest av- 
eraged intensity over the search area for simultaneously all the 
wavelengths. 
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Fig. 4. Schematic illustration of the wavelengths and bandwidth 
of study. Ai and Ao represent the limit wavelengths of the spectral 
bandwidth. 



3.4. Computation of the broadband image in plane D 

Once the parameters optimized, we can calculate the broad- 
band image obtained with the CAPM coronagraph using Eq.©. 
Since o>\j, u)2,a and fix are estimated for five wavelengths, the 
monochromatic images at those wavelengths can be computed. 
To sum more than five frames for our final broadband image, 
we need to know the parameter values at other wavelengths. An 
interpolation of the values is realized for each apodization pa- 
rameter. After several tests, we decide to use a second degree 
polynomial function to fit the parameter values, more represen- 
tative of the behavior of actual materials than higher order func- 
tions. In the following, all our interpolations are performed with 
a second degree polynomial function. Once this defined, we can 
estimate the values of the parameter at a given wavelength and 
compute its frame. We then sum several monochromatic images 
with small spectral separation between two successive frames 
and obtain our broadband image in plane D. 
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4. Performance over large spectral bandwidths 

In this section, we study the behavior of the CAPM corona- 
graph and its theoretical performance for different spectral band- 
widths. This investigation is made for a clear circular aperture 
but it can be generalized to any kind of entrance pupil geom- 
etry. The current study will be particularly interesting in the 
context of some exoplanet imagers under consideration, such 
as 1-2 m class off-axis spac e telescopes equipped with an in- 
ternal coronagraph (see e .g., iGalicher et al.ll2Q lot iTrauger et al.l 
120 lOt iGuvon et al.ll2010ch . In addition, the results reached with 
our concept are comp ared to those previously achieved by 
ISoummer et al. (2003b) with the DZPM coronagraph. This al- 
lows us to underline the contrast gain obtained with the use of a 
colored apodization instead of a grey one. 

4.1. Colored apodization 

Let us first observe the shape of the apodization reached for the 
CAPM coronagraph at different wavelengths. Figure [5]panel (a) 
shows the radial intensity profiles of the colored apodization at 
five wavelengths for a 25% bandwidth (H-band). We can notice 
the similarity between all the displayed profiles. To disentan- 
gle them, the difference between the transmission functions at 
a given wavelength A and at Aq has been drawn on the panel 
(b), emphasizing small difference in absolute value existing be- 
tween different profiles. The intensity difference is not larger 
than 3% of the normalized intensity peak but it reveals the chro- 
matic aspect of our apodization used for the CAPM corona- 
graph. The apodization throughput evolves slightly from 57.9% 
(A= 1.485 ^im) to 55.3% 1.815/mi). 
The panels (c) and (d) of Figure [5] show the phase ip of the phase 
apodization O,, = e lip in a similar way. The phase difference be- 
tween the phases at a given wavelength A and at Ao is not larger 
than 0. 1 rad but it emphasizes the chromatic aspect of the phase 
apodization <t>„. in the CAPM coronagraph. 

4.2. Residual intensity 

Much information can be extracted from the monochromatic in- 
tensity profiles in the Lyot and image planes, see Figure[6] These 
curves have been obtained considering the optimized apodizer 
showed in Figure[5] The minimization of the residual broadband 
intensity inside the search area (2-10/lo/D range) has driven the 
optimization of our CAPM concept, leading to coronagraphic 
images with high intensity peak and low sidelobe level, see 
Figure [6] bottom plot. In parallel, the energy distribution inside 
the Lyot stop exhibits apodized intensity profiles, see Figure [6] 
top plot. Observing all the wavelengths, we note that the bright- 
ness of the sidelobes beyond 6Aq/D angular separation follows 
the intensity value found at the pupil edge. For sidelobes in the 
2-6 Aq/D range, this assertion is no longer valid because of the 
combination of several effects, including oscillations of the am- 
plitude for some wavelengths in the Lyot plane as well as the 
presence of a phase apodization in the entrance pupil. 

4.3. Optimized parameters 

Following our optimization criteria, we estimate the different pa- 
rameters of the CAPM coronagraph for different spectral band- 
widths. The obtained values are given in TableQ] In addition, the 
colored apodization parameters oji j, a>2 t A and /3a are plotted as a 
function of the wavelength for different bandwidths in Figure|7] 
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Fig. 5. (a): radial profile of the colored intensity apodization for 
the CAPM coronagraph at different wavelengths A. (b): radial 
profiles of the difference between the intensity apodization at a 
given wavelength A and that obtained at the central wavelength 
Aq (here 1.650 fj.m). (c) and (d) reproduce the panels (a) and (b) 
respectively, for the phase given by the phase apodization <t> n ,. 
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Fig. 6. Radial profiles of the residual intensity with the CAPM 
coronagraph, within the re-imaged pupil in plane C (top) and 
in the final image plane D (bottom) at different wavelengths A. 
Profiles have been normalized at the origin to directly read the 
coronagraphic extinction ratio at a given angular distance from 
the on-axis star image. 



4.4. Contrast results and discussion 

In Figure [8] we draw the radial intensity profile of the im- 
ages achieved with the CAPM coronagraph for different spec- 
tral bandwidths. The plotted curves are extracted from polychro- 
matic images computed with N= 1 + 100 AA/Aq monochromatic 
frames with a separation corresponding to 1% band. The val- 
ues of the averaged intensities reached for each band and at sev- 
eral angular separations are given in Table Q] Averaged inten- 
sities of 2.2 ■ 10~ 7 and 2.2 • 10" 8 are theoretically reached at 
3 Aq/D and 5 Aq/D respectively, with our coronagraphic device 
for a 25% bandwidth. As we can notice in Figure|7] the interpo- 
lation function is not optimal for the largest bands, leading to a 
loss of coronagraph performance when the optimized parameters 
are replaced by the interpolated values: the averaged intensity at 
5 A Q /D degrades from 3.8 ■ 10~ 8 and 7.1 • 10~ 8 to 5.0 ■ 10" 8 and 
1 .3 • 10~ 7 for the 40% and 50% bands respectively. 
To go further in the analysis, we also determine the averaged 
intensity at 5 Aq/D achieved by the coronagraph at each wave- 
length of a given spectral bandwidth case, see Figure [9] Each 
curve is plotted with the same 1 % sampling as described above. 
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Fig. 7. Optimized values of the amplitude apodization parameter 
(jD\ j (top plot), (l>2,a (middle) and the phase apodization param- 
eter /3a (bottom) as a function of the wavelength for different 
bandwidths. A second-degree polynomial interpolation of the 
points is represented for each parameter and bandwidth. 



It can be noticed from this plot that a 2 • 10 averaged in- 
tensity at ~ Aq can theoretically be achieved in the case of a 
25% bandwidth. Another interesting aspect about reachable av- 
eraged intensities concerns the spectral range. As the bandwidth 
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decreases below 20%, a single hole in the intensity curve ap- 
pears around Aq, surrounded by increasingly steep edges. In the 
monochromatic case, the coronagraph becomes identical to the 
ARPM coronagraph as we could expect. 

Finally in Figure [TO] we compare the results achieved with 
the CA PM coronagraph with those obtained by Sou mmer et alj 
(2003b) with the DZPM coronagraph. It can be noticed that a 
~2.5 mag contrast gain at 3Ao/D is reached with our device 
compared with the first DZPM coronagraph design for a 20% 
and a 40% bandwidths. Similar improvements of the reachable 
averaged intensities have been observed at other angular sepa- 
rations (5 and 7 Aq/D), confirming the substantial gain obtained 
for the CAPM coronagraph by replacing the grey by a colored 
apodization. 
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achieved with the CAPM coronagraph for different sizes of the 
spectral bandwidth. 
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Fig. 9. Theoretical averaged intensity at 5 Aq/D obtained with 
the CAPM coronagraph as a function of the wavelength for dif- 
ferent bandwidths. 
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gular separations from the main optical axis with the DZPM and 
CAPM coronagraphs as a function of the bandwidth. 



5. Design proposal for the colored apodization 

The realization of a colored apodizer is a key point for the im- 
plementation of the CAPM coronagraph. This requires careful 
control of the transmission function of the apodizer both spa- 
tially across the pupil and spectrally over a wide range of wave- 
lengths. We expose here some of the options currently investi- 
gated for the fabrication of this chromatic apodization. 
One option consists of using an absorbing material with an ap- 
propriate chromatic variation of the absorption coefficient, for 
instance a colored-glass filter. The required transmission func- 
tion for the apodizer is obtained at all the wavelengths by care- 
ful thickness control of the absorbing material. Appropriate ma- 
terial choice would lead to a negative lens of up to a few mm 
edge thickness. The optical power of this lens could possibly be 
included in the optical design of the instrument, but in a more 
general option it would need to be compensated by a positively 
powered non-absorbing lens. 

Another approach for the colored apodizer that we have consid- 
ered is to use a thin film of a highly absorbing material such as 
metal, deposited with an appropriate radial thickness profile onto 
a silica substrate. The chromaticity of the apodizer is now related 
to the chromatic variation of the complex refractive index of the 
material. 

Solutions based on metallic thin film, alloy layer or else High- 
Energy Beam Sensitive (HEBS) glass have already been inves- 
tig ate d to produce achromatic band-l i mited masks in the vis- 
ible (Sidickl 120071 iMoodv et all 12008; Balasubjamanian|[2008) 
or manufacture achr omatic apodizers of t he AP L C coronagraph 



for V ET-SPHERE dGuerri et alJ 120081. 1201 ll 
201 ll) and Gemini Planet Imager (ISoummer et al 



ICarbillet et all 

200' 



Sivaramakrishnan et al . 2009). Our study is focused here on the 
search of metals to design chromatic components instead of grey 
masks. 

In the following, we detail this approach for the realization of 
the colored apodization. The performance of the CAPM corona- 
graph is estimated with this physical design and then compared 
to the results with the theoretical model of colored apodizer de- 
scribed above. We assume a clear circular aperture and a 25% 
bandwidth centered at Aq = 1 .65 fim (H-band) for our optimiza- 
tion. 
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Table 1. Values of the CAPM coronagraph parameters optimized for different spectral bandwidths and averaged intensity results at 
different separations. The DZPM phase steps are given in optical path difference (OPD). Five wavelengths from A\ to /I5 have been 
chosen within the spectral bandwidth AA, centered on A?, = Aq, and equally spaced from one to another by AA/5. 



Parameters 






Bandwidth in % 








5 


10 


20 


25 


30 


40 


50 


d\ in Ao/D 


0.877 


0.881 


0.879 


0.874 


0.880 


0.883 


0.876 


d 2 in Ao/D 


1.424 


1.430 


1.438 


1.445 


1.455 


1.479 


1.497 


OPD! in A 


0.311 


0.312 


0.311 


0.309 


0.312 


0.313 


0.314 


OPD, in A Q 


0.675 


0.677 


0.679 


0.678 


0.684 


0.690 


0.691 


oji j at Ai 


-2.401 


-2.390 


-2.422 


-2.420 


-2.402 


-2.332 


-2.182 


wi,^ at A2 


-2.402 


-2.391 


-2.371 


-2.368 


-2.374 


-2.404 


-2.417 


a) 1A at /1 3 


-2.424 


-2.424 


-2.420 


-2.404 


-2.397 


-2.392 


-2.320 


Oil A at ^4 


-2.471 


-2.488 


-2.510 


-2.509 


-2.522 


-2.502 


-2.417 


(jj i A at A 5 


-2.484 


-2.536 


-2.628 


-2.668 


-2.674 


-2.741 


-2.705 


u) 2 j at li 


2.289 


2.282 


2.474 


2.508 


2.621 


2.823 


2.924 


<jJ2,a at Ai 


2.275 


2.246 


2.170 


2.135 


2.207 


2.349 


2.454 


cu 2 j at A 3 


2.340 


2.328 


2.282 


2.193 


2.176 


2.118 


1.847 


cu 2A at A 4 


2.492 


2.534 


2.567 


2.538 


2.582 


2.493 


2.225 


U) 2A at A 5 


2.530 


2.679 


2.954 


3.072 


3.104 


3.305 


3.241 


Pa at Ai 


-0.229 


-0.235 


-0.245 


-0.251 


-0.259 


-0.272 


-0.265 


Pa at A 2 


-0.229 


-0.231 


-0.233 


-0.234 


-0.240 


-0.251 


-0.255 


Pa at A 3 


-0.228 


-0.230 


-0.227 


-0.224 


-0.225 


-0.224 


-0.215 


Pa at A 4 


-0.229 


-0.230 


-0.228 


-0.226 


-0.226 


-0.221 


-0.214 


Pa at A 5 


-0.229 


-0.232 


-0.236 


-0.239 


-0.240 


-0.249 


-0.252 


level at 1.5 Ao/D 


8.8- 10- 5 


9.5- 10- 5 


1.110 4 


1.3-10 4 


1.4-10" 4 


1.7-10~ 4 


2.0- 10~ 4 


level at 2.0 /V-D 


6.6- 10" 6 


7.4- 10~ 6 


9.9- 10" 6 


1.210 5 


1.310 5 


1.8-10" 5 


2.4- 10~ 5 


level at 3.0 Ao/D 


5.2-10~ 8 


7.6- 10~ 8 


1.610" 7 


2.2-10" 7 


2.4- 10~ 7 


3.8-10- 7 


7.2- 10~ 7 


level at 5. (M AD 


1.110 8 


1.2- 10-* 


1.9-10" 8 


2.2-10~ 8 


2.9- 10~ 8 


3.8-10" 8 


7.1-10 8 


level at 7.0 A /D 


5.9-10- 9 


6.3-1Q- 9 


8.3-10" 9 


9.0- 10-" 


1.2-10" 8 


1.410 8 


2.5- 10~ 8 



5. 1. Assembly based on a metallic thin film on a silica 
substrate 

For the manufacturing of our colored amplitude apodization O a , 
we investigate an assembly based on a metallic thin film de- 
posited on a silica substrate with complex refractive index N m 
and N s respectively, see Figure [TT] At the boundaries of our as- 
sembly, we consider an absorption-free medium of refractive in- 
dex A^o- The assembly is assumed to work in transmission at nor- 
mal incidence. 

The substrate is supposed to be an absorption-free medium, thick 
enough to neglect the internal reflections. Referring to the metal- 
lic thin film, the refractive index N m presents a complex part, 
called extinction coefficient k, which is related to the absorption 
a m within the metal layer by a m - Ank/A. Increasing the metal 
thickness leads to a combined enhancement of the absorption 
and reflection of the assembly and therefore, a reduction of its 
transmission. A careful control of the metallic thin film thick- 
ness allows us to adjust the amplitude transmission function of 
our assembly. The chromaticity of the amplitude apodization is 
mainly due to the dependence on wavelength of the metal ex- 
tinction coefficient. 

For the present concept, the amplitude transmission function of 
our assembly <!>/, is defined as the product of the transmittances 
through the metallic thin film and the substrate r m and r s respec- 
tively. The amplitude apodization can therefore be written as: 

*&(r,A)=r m (r,A)r,(A), (12) 

For the sake of clarity, we omit the dependence on wavelength of 
the refractive indices and the parameters B, C and 5,„ which ap- 
pear in the following equations. The transmittances in amplitude 



r m and t s can be expressed as: 

T m (r,A)=2N KN B + C), 

t s (A)=2N s /(N s +N ), 

where B and C can be identified as the normalized electric 
and magnetic fie ld amplitudes at the front interface respectively 
(lMacLeodlll99ll) . In our case, B and C can be calculated with 
the following equation: 

B\_l cos6 m (ism8 m )/N m \( 1 \ 
CJ \iN m sm6 m cosc5 ffl J\N S J ^ ' 

with 6,„(r, A) = 2iiN m t{r)l A the phase shift experienced by the 
wave as it traverses a distance t normal to the substrate boundary. 
The complex apodization <t> of the CAPM coronagraph with the 
physical design writes as: 

<D = <D fe • <b w . (15) 

The metallic thin film also induces a radially variable phase de- 
lay, leading to complex values evolving with r for the transmis- 
sion function 0^. Since our concept of CAPM coronagraph also 
contains a phase apodization <D, V , the total phase of the complex 
apodization <t> in the entrance pupil plane A will be the combi- 
nation of the phase terms of <$> w and O/,. 

5.2. Choice of the metal 

Preliminary studies have been realized to choose the metal for 
the physical model, see Figure[T2] Five metals (gold, silver, alu- 
minum, nickel and chromium) have been considered for this pur- 
pose while five radial distances have been investigated to deter- 
mine the optimal thickness of the metallic thin film at different 
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Fig. 11. Schematic illustration of the assembly, composed of a 
metallic thin film deposited on a substrate. The scheme is not at 
scale. 



positions in the pupil. For each metal and radial distance, the 
thickness of the metallic thin film has been estimated such that 
the assembly reaches the same intensity transmission as the the- 
oretical model at Aq observed in Figure [5] B ulk re fractive in- 
dices for the metal layer, obtained from Palik| dl985h . have been 
used for the implementation of our simulations. The thicknesses 
reached for each metal and radial distance are given in Table|2] 
Using silver or nickel provides the best fit to the theoretical 
model at small or large radial distances, respectively. However, 
even with those metals, the theoretical profile of the colored 
apodization is not reached at all the wavelengths of study. In 
the following, we re-optimize the CAPM coronagraph, using a 
nickel apodizer. 

With reference to the apodizer model of Eq. ©, the thickness 
profile t of the metal is optimized using the following expres- 
sion: 

[ Is / " 4> 



t(r) = t c 



r 4 - 



(16) 



where t c = 1 nm is the reference thin film thickness and x\ 
and xi are the parameters related to the terms in r 2 and r 4 
respectively of the thickness profile. This assembly provides a 
maximum amplitude transmission at the center of the pupil or at 
the edge of the central obscuration for rj + 0. 



Table 2. Thickness of the different metals for the assembly at 
different radial distances used in Figure [T2l 



Radial distance 


Metallic thin film thickness in nm 


Au 


Ag 


Al 


Ni 


Cr 


r = 0.125 


1.044 


1.237 


0.324 


0.485 


0.557 


r = 0.250 


2.822 


3.228 


0.977 


2.070 


2.449 


r = 0.375 


5.081 


5.733 


1.844 


4.952 


6.092 


r = 0.500 


8.177 


9.148 


3.062 


9.666 


12.403 



5.3. Results 

The optimized parameters for the phase mask and apodizer are 
reported in Table [3] while the resulting thickness profile t of the 
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Fig. 12. Transmission values in intensity obtained for the theo- 
retical model and the assembly with different metals as a func- 
tion of the wavelength. Values are provided for different normal 
radial distances. Thickness for each material and radial distance 
are given in Table [2] 



nickel film and the parameter /3a associated to the phase apodiza- 
tion are plotted in Figure [13] Regarding t displayed in the top 
plot, a maximum thickness of 9. 1 nm is reached at the pupil edge 
for the metal layer, proving to be consistent with the 9.7 nm value 
observed at r — 0.5 in the tests for the choice of the metal and 
reported in Table [2] Concerning the /3 a curve in the bottom plot, 
we can note some discrepancy between the values obtained for 
the CAPM coronagraphs working with the theoretical model and 
the physical design. This difference can be explained by the pres- 
ence of the additional phase shifting introduced by the metallic 
thin film of the assembly, leading to a re-adjustment of the pa- 
rameter /3a of the chromatic defocus during the optimization of 
the CAPM coronagraph with the physical design. 
The radial intensity profiles of the apodizer for a 25% bandwidth 
(H-band) are represented at five wavelengths in Figure [T4lpanel 
(a). While the profiles are quite similar to the profiles shown in 
Figure [5] the difference between them, shown in panel (b), are 
smaller than for the optimal profiles. The apodization through- 
put of the assembly ranges from 59.5% (A — 1 .485 ^m) to 58.8% 
(A = 1.815 //m), which is slightly higher than the theoretical col- 
ored apodizer. The phase profiles, shown in panel (c), are some- 
what different from those of the optimal apodizer in Figure [5] 
indicating the presence of a fourth-order term (spherical aberra- 
tion) in addition to the second order term. The differential phase 
profiles in panel (d) are remarcably similar to those of the opti- 
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mal apodizer, both in shape and amplitude. 
In Figure Q3J we represent the intensity profiles obtained with 
the CAPM coronagraph in the presence of the theoretical col- 
ored apodizer and our physical design. The profiles are almost 
identical, underlining the good relevance of using a nickel layer 
to realize our colored amplitude apodization. The intensity lev- 
els estimated at different separations from the star are thus very 
close for the two models of colored amplitude apodization, see 
Table |U For instance, the intensity levels at 5 Aq/D are 2.6 ■ 1CT 8 
and 2.2 • 1CT 8 , respectively. An optimized design using a nickel 
thin film therefore allows us to reproduce the contrast results 
obtained for the optimal CAPM coronagraph. The nickel layer 
presents a maximum thickness of 9.1 nm with an accurate con- 
trol at atomic level. To manufacture such a design, atomic layer 
deposition (ALD) represents a very attractive technique, widely 
used in microelectronics, since it enables the production of very 
thin, conformal film with thickness control and composition 
of the films at the a tomic level for a large range of materials 
dLeskel & Ritalall2002l) . 

Further studies and tests will be required for the manufactur- 
ing of such chromatic apodizers. Finally, options based on more 
elaborate, radially variable thin film stacks, possibly containing 
both dielectric and metallic materials, are also being considered. 



Table 3. Values of the parameters for the CAPM coronagraph 
optimized with an assembly using nickel thin film for the col- 
ored apodization. The DZPM phase steps are given in optical 
path difference (OPD). Five wavelengths from A \ to A5 have been 
chosen within the spectral bandwidth AA, centered on A3 = Ao, 
and equally spaced from one to another by AA/5. 



Parameters 


25% bandwidth 


d] in A /D 


0.868 


d 2 in A /D 


1.478 


OPDi in A 


0.338 


OPD 2 in A Q 


0.716 


Xi 


28.9 


Xi 


30.2 


Pa at -l. 


-0.320 


Pa at A 2 


-0.318 


Pa at A 3 


-0.319 


Pa at A 4 


-0.331 


Pa at ^5 


-0.358 



Table 4. Intensity values provided by the CAPM coronagraph 
with the theoretical model and the physical design for the col- 
ored apodization. Results are given for a 25% bandwidth cen- 
tered at Aq = 1.650/mi (H-band). 



level at 
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1.4- 10- 7 
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Fig. 13. Top: thickness profile of the nickel layer as a function of 
the wavelength. Bottom: optimized values of the phase apodiza- 
tion parameter /3a for the CAPM coronagraph in the presence of 
the theoretical model and the physical design as a colored am- 
plitude apodization. 



5.4. Implementation of the chromatic defocus 

The colored apodizer also includes phase apodization <!>„ with a 
chromatic parameter [3 a- While grey phase apodization can sim- 
ply be realized by defocussing the mask, we address here the 
possibility to produce a chromatic defocus. Our numerical opti- 
mization shows that parameter^ varies as A 2 for the theoretical 
model and our physical design, see Figure[l3] A possible way to 
realize a 2 nd degree polynomial chromatic defocus consists in us- 
ing a powerless doublet at the entrance pupil plane. This doublet 
will present the following characteristics: infinite focal length at 
/to, a null 1" order chromatic dispersion and a constant and non 
null 2 nd order chromatic dispersion, also called anomalous dis- 
persion. This can be achieved by combining glasses which have 
the same refractive indices n\ and «2, the same Abbe numbers 
v\ and V 2, but different partial dispersions. This work follows 
iDuplovl (120061) who has recently realized studies of apochro- 
matic lens design by considering the partial dispersion proper- 
ties of materials from standard glass catalogs rather than special 
anomalous-dispersion glasses. We have performed preliminary 
tests of this approach in the visible range using ZEMAX® ray 
tracing code and tabulated glass data. Encouraging results for 
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Fig. 14. (a): radial intensity profile of the colored apodizer, made 
with an assembly using nickel layer, for the CAPM coronagraph 
at different wavelengths A. (b): radial profiles of the difference 
between the intensity apodization at a given wavelength A and 
that obtained at the central wavelength Aq (here 1.650jum). (c) 
and (d) reproduce (a) and (b) respectively for the phase given by 
the phase apodization O w . 



the realization of our phase apodization were found with a com- 
bination of Schott glass materials N-FK9 and K5, see Figure[T6l 
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Fig. 15. Radial intensity profile of the coronagraphic images 
achieved with the CAPM coronagraph using the theoretical 
model and the physical design for the colored apodizer. 



Indeed, chromatic behavior in A 2 is obtained with this doublet, 
corresponding to the chromatic defocus expected for the phase 
apodization of our CAPM coronagraph. Combination of other 
glass materials are currently under study to achieve adequate 
chromatic defocus in other spectral bands. 



6. Sensitivity analysis 

After optimizing the CAPM coronagraph for achromatization, 
we investigate here several effects that may alter the performance 
provided by this coronagraphic system. In the following sensi- 
tivity analysis, we consider a CAPM coronagraph optimized for 
a 25% bandwidth (H-band) observing a point-like source, ex- 
cept for the star angular size study, in the case of a clear circular 
aperture. 



6.1. Performance criteria 

Contrast performance is judged by the averaged intensity esti- 
mated at 5 Ao/D from the main optical axis in the image plane 
and for an annulus of width Ao/D. 

To quantify the capacity to work close to the star, the Inner 
Working Angle (IWA) is often used, defining the minimum an- 
gular distance at which the throughput of a companion is larger 
than 50%. Phase mask coronagraphs offe r very small IWA, typi - 
cally ~ Ao/D, as recalled for instance bv lMartinez et alJ (|2008). 
This is the case for the CAPM coronagraph as illustrated in 
Figure [T7] Unfortunately, IWA is not related to the notion of 
contrast and therefore, the intensity of close-in companions de- 
tectable at a given angular distance from the star cannot be de- 
termined with this benchmark. 

We propose a new criterion, named Inner Detectability Angle 
(IDA), which represents the compromise between IWA and 
achievable contrast. IDA is defined as the minimum angular dis- 
tance a from the main optical axis for which the ratio between 
the normalized intensity of our coronagraphic star image and the 
planet transmission flux T p is less than a given value, here set to 
1(T 6 : 

I D (a) 1 



T P (a) 



< 10" 



(17) 
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Fig. 16. Top: design of a lens doublet combining Schott glass 
materials N-FK9 and K5. Bottom: focal shift versus wavelength. 
This result has been obtained assuming a longitudinal focal shift 
reached in the visible range with our values of /3a- A 2 nd order 
polynomial chromatic focal shift with a 5 mm amplitude range 
has been achieved. 



where 1$ denotes the intensity peak value in the absence of coro- 
nagraph mask. In Figure [TT] we plot the planet transmission 
profile as a function of the planet position in the presence of a 
CAPM coronagraph optimized for a 25% bandwidth. From this 
curve, we deduce that the planet transmission can reasonably be 
taken equal to 1 for distances larger than 1.5 Aq/D. IDA is ex- 
pected to be larger than IWA, but estimating IDA allows us to 
know at which angular distance companions, 10~ 6 fainter than 
their host star, could be detected with our CAPM coronagraph. 
We study the evolution of IDA for the different cases described 
below. 



6.2. Effects of aberrations 

Wavefront errors constitute a major issue for stellar corona- 
graphs as they limit the performance provided by these diffrac- 
tion suppression systems. This crucial point of aberrations then 
has to be addressed in the framework of a coronagraphic planet 
imager. In the case of a ground-based telescope with an XAO 
system, as in the case of a well-corrected space instrument, the 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 
Planet position from the star in A /D unit 

Fig. 17. Planet transmission profile achieved with the CAPM 
coronagraph as a function of the planet position. 



residual quasi-static instrumental aberrations are known to affect 
the intensity level reached by a downstream coronagraphic de- 
vice. A performance estimate of the CAPM coronagraph in the 
presence of aberrations is therefore required to know the sensi- 
tivity of our device. In the following, we describe the behavior 
of our coronagraphic concept in the presence of low-order aber- 
rations (tip-tilt, defocus, astigmatism, coma, and spherical aber- 
ration) since they represent the main contribution in terms of 
wavefront errors. In addition, we analyze the impact of random 
aberrations caused by real optical elements on the CAPM coro- 
nagraph performance. This aims to determine the limits of our 
system in the presence of light scattered by the surface rough- 
ness of optical elements. 

The low-order aberrations are decomposed here following the 
Zernike polynomials basis. Their amount is expressed in Zernike 
rms coefficients and given in Ao unit. The random aberrations 
considered here are related to the wavefront errors introduced by 
real optical elements. Good quality optical surfaces are found to 
exhibit a power spectral density (PSD) function proportional to 
where v de notes the spatial frequency i n the pupil, as recalled 



,,-2 



for instance bv lRoddier & Roddierl (Il997l) . To realize our study 
with random aberrations, we introduce numerically a randomly 
aberrated phase screen following this PSD, in the entrance pupil 
plane A of our coronagraph. 

6.2.1. Qualitative analysis 

Figure [18] gives a set of the images theoretically obtained with 
the CAPM coronagraph. These frames are represented for differ- 
ent values of rms wavefront errors and each aberration consid- 
ered here. The shape of the coronagraphic image modifies and 
evolves quickly as the amount of the rms wavefront errors in- 
creases. 

In Figure [19] we draw the radial intensity profiles of the corona- 
graphic images obtained with the CAPM coronagraph for differ- 
ent rms wavefront error values of defocus (top plot), spherical 
aberration (middle plot) and random aberration (bottom plot). 
Through the evolution of the profile with the amount of a given 
aberration, one can appreciate the variation of the IDA and the 
intensity level reached with the CAPM coronagraph. A dashed 
line representing the 10~ 6 intensity level set for IDA is also dis- 
played on the plot. 
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Fig. 18. Set of coronagraphic images achieved with the CAPM coronagraph and computed numerically for different aberrations. All 
the images are normalized to the peak of the non coronagraphic image. The pixel scale is the same for all the images and the field 
of view is represented with a linear range of 22 Aq/D, 



We note that while both low-order and random aberrations af- 
fect all angular separations, the low-order profiles remain paral- 
lel with the non-aberrated profile, following a a~ 3 law, as in the 
case of the classical Airy pattern. In the random aberrations case, 
on the other hand, the profile changes towards a or 2 profile, cor- 
responding to the PSD profile defined above. We conclude that 
mechanism leading to performance degradation is different in 
the two cases. While edge diffraction dominates in the low-order 
case, wavefron t diffraction, as modeled by a sum of phase grat- 
ings (iGoodmanll 19961) dominates in the random case. 



6.2.2. Intensity & IDA 

In Figure [20] top plot, we represent the averaged intensity at 
5 Aq/D achieved by the CAPM coronagraph as a function of 
the rms wavefront errors value. The averaged intensity goes be- 
yond a 10~ 7 value for rms wavefront errors value larger than 
~ 3 • 10" 3 Aq and ~ 9 • 10 Aq for low-order and random aber- 
rations respectively. At greatest aberration values, the reachable 
averaged intensity increases approximately as a 2 nd degree poly- 
nomial function for all the aberrations. 

Let us note that in the case of a spherical aberration, a small and 
local improvement of the coronagraph performance is observed 
for a wavefront error of ~ 7 • 10~ 4 /lo rms with respect to the 
reference averaged intensity value in the absence of aberration. 
We deduce that an additional phase apodization in the form of 
spherical aberration can slightly improve the contrast gain pro- 
vided by this device but the investigation of this effect is not 
within the scope of this paper. However, the gain at 5 Ao/D im- 
plies a non desirable increase in intensity of the first Airy ring, 
see the 10" 3 ^o 

rms profile in Figure [19] middle plot. 
To complete the analysis, we represent the IDA theoretically 
achieved with the CAPM coronagraph as a function of the rms 



wavefront errors value, see Figure [20] bottom plot. IDA remains 
below 3 Ao/D for rms wavefront errors lower than 3 • 10" 3 /lo 
and 2 • 10~ 3 Aq for the low-order and random aberrations respec- 
tively. At larger values, the IDA evolves approximately as a 2/3 
and a second degree polynomial function for the low-order and 
random aberrations respectively. 



6.2.3. Discussion 



The intensity level and IDA provided by CAPM coronagraph 
are all the more altered as the aberration order is higher. 
These results show that at a given value of rms aberration, the 
sensitivity of our coronagraph increases with the order of the 
aberration. They also underline the necessity of working with 
very high quality optical elements, typically better than 10~ 3 Ao 
rms wavefront error, to enable high performance of the CAPM 
coronagraph for a very high contrast imaging as well as a very 
small IDA. 

Several studies of sensitivity to low-order aberrations 
have been performed for d i fferent kinds of coron agraph 
(e.g IShaklan & Greenl l2005t [sl varamak rishnan et aTl [2005; 
Belik ov et alJ l2006t lMawetetalJl2010h . No general behavior 
of the coronagraphs to low-order aberrations emerge from 
the previous analysis since each concept presents a specific 
response to these wavefront errors. In addition and to the best 
of our knowledge, no advanced study has been reported on the 
response of circular phase mask coronagraphs to low-order 
aberrations. 
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Fig. 19. Radial intensity profiles of the coronagraphic images 
theoretically achieved with the CAPM coronagraph for different 
rms wavefront errors values of defocus (top plot), spherical aber- 
ration (middle) and random aberration with PSD defined above 
(bottom). 
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Fig. 20. Theoretical averaged intensity at 5 Aq/D (top plot) and 
IDA (bottom plot) achieved with the CAPM coronagraph for dif- 
ferent aberrations. 



6.3. Effects of mask imperfections 
6.3.1. Mask roughness 

Effects of the mask roughness have to be analyzed since they 
could alter the performance of the CAPM coronagraph. Many 
glass substrates with excellent flatness (better than A/20, see for 
instance Newport ©10QW40-30 glass) are available and useful 
for our DZPM manufacturing. To lead our study, we therefore 
assume a perfectly flat substrate of the DZPM except where the 
phase disk and annulus are machined into it. So, only effects of 
surface roughness present inside the mask area of diameter d.2 
have been considered. 

In Figure [21] top plot, we draw the radial profile of the corona- 
graphic image for different values of rms mask roughness. All 
the coronagraphic profiles follow a or -3 law as that of the Airy 
diffraction pattern does and their intensity level increases with 
the amount of rms mask roughness. 

Following this result, we also represent the averaged intensity at 
5 Ao/D reached by the CAPM coronagraph as a function of the 
rms mask roughness, see Figure|22]top plot. It can be noticed that 
a 10 intensity level is provided with the CAPM coronagraph 
for rms mask roughness values lower than ~ 0.12/1(>. Beyond 
this value, the intensity increases and follows the behavior of 
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a 2 nd degree polynomial function until it reaches a plateau at 
2 TO 5 , corresponding to the non-coronagraphic case, for rough- 
ness greater than 1.00 Aq rms. 

In addition, the IDA profile is displayed in Figurel22lbottom plot. 
A 3 Aq/D IDA is preserved for rms mask roughness value lower 
than 0.15 /in. For larger roughness values, the IDA increases 
and follows a 1 " degree polynomial function, again reaching a 
plateau at 1.00 /In rms. 

For 0.12 Aq rms mask roughness, the corresponding values are 
198.0 nm and 76.0 nm at A = 1650 nm (H-band) and at 
/In = 633 nm (visible) respectively. As a reference, the rough- 
ness obtained with ou r second-generation RRPM is 0.8 nm rms 
dN'Diave et al.ll20Toh . Since the same accuracy can be obtained 
for the manufacturing of DZPM as for that of the RRPM, we can 
deduce that mask roughness will not constitute a major issue for 
the CAPM coronagraph. 



6.3.2. Mask transition zone width 

The DZPM is theoretically a hard-edge mask. However in prac- 
tice, some transition zone widths are found in a manufac- 
tured prototype as noticed with the second-generation RRPM 
dN'Diave etalJl20Toh . We propose to study the impact of the 
DZPM transition width on the CAPM coronagraph performance 
by replacing vertical transitions by slopes in the mask shape. The 
same width is assumed here for the inner and outer parts of the 
DZPM and expressed in mask outer part diameter &i. The tran- 
sition zone is illustrated in Figure [2] 

Radial profiles of the coronagraphic images are displayed for 
different mask transition widths in Figure[2T|bottom plot. All the 
profiles evolve as a 3 rd degree polynomial function like the Airy 
diffraction pattern. As expected, the intensity values of the coro- 
nagraphic image increase as the DZPM transition zone width 
becomes larger. 

In Figure[22]top plot, we display the averaged intensity at 5 Aq/D 
provided by the CAPM coronagraph as a function of the DZPM 
transition zone width. An intensity level of 10~ 7 is preserved for 
a mask transition width smaller than 0.014^2- At larger widths, 
the averaged intensity increases but slower than a 2 nd degree 
polynomial function. 

Figure [22] bottom plot shows the evolution of IDA achieved 
with the CAPM coronagraph as a function of the mask transi- 
tion width. We can notice that IDA remains below 3 Aq/D for 
a mask transition width lower than 0.028 d^. Beyond this value, 
IDA evolves as a 1/2 degree polynomial function, showing the 
low sensitivity of this criteria to an increase of the mask transi- 
tion width. 

As a summary, a 1% mask transition width allows us to pre- 
serve a 10~ 7 intensity level and good IDA with the CAPM coro- 
nagraph. This value is similar to that measured for the second- 
generation RRPM wh ich was less than 1 % of the mask diameter 
dN'Diave et alJ 2010) and below 1 //m. The same accuracy is ex- 
pected concerning the manufacturing of the DZPM. A 1 % mask 
transition width corresponds physically to 1 jum for a 100 fim 
DZPM diameter. Hence, the focal ratio F of the coronagraphic 
layout will range as follows: F ~ 61 at A — 1.65 /mi, F ~ 100 at 
A = 1 jt/m and F ~ 200 at A = 0.5 //m. Below a 1 fim wavelength, 
F constitutes a severe constrain on the design of a planet im- 
ager working in the visible and in which a CAPM coronagraph 
could be inserted. To overcome focal ratio issues in the visible 
spectral range, manufacturing DZPM with a diameter and mask 
transition width smaller than lOO/im and 1 fim respectively is 
required. 
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Fig. 21. Radial intensity profiles of the coronagraphic images 
theoretically achieved with the CAPM coronagraph for differ- 
ent rms values of mask roughness (top plot) and for different 
sizes of mask transition width, expressed in DZPM diameter dt 
(bottom). 



6.4. Lyot stop size and misalignment 
6.4.1. Lyot stop size 

Pupil undersizing in the Lyot plane is often required in the con- 
text of instrument design to allow alignment and manufacturing 
tolerance. The sensitivity of the CAPM coronagraph to this pa- 
rameter has therefore to be taken into account. We analyze the 
impact of a Lyot Stop (LS) size reduction on the contrast pro- 
vided by the CAPM coronagraph. Two CAPM coronagraphs op- 
timized with a 90% and a 100% LS are considered here. 
In Figure [23] the averaged intensity at 5 Aq/D achieved by the 
CAPM coronagraph is represented as a function of the pupil size 
for both designs. In each case, the best contrast provided by the 
CAPM coronagraph is almost obtained for the LS size for which 
the coronagraphic device has been optimized. However, the be- 
havior of the CAPM coronagraphs differs from each other. In the 
case of the CAPM coronagraph optimized with a 100% LS, the 
intensity level increases from 2 TO -8 to 2 TO -7 as we reduce the 
size of the LS from 1 D to 0.85 D. On the opposite, in the case of 
the CAPM coronagraph optimized with a 90% LS, the intensity 
level varies slowly from 3 TO" 8 to 5 T0~ 8 for LS sizes diverg- 
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Fig. 22. Theoretical averaged intensity (top plot) and IDA (bot- 
tom plot) achieved with the CAPM coronagraph as a function of 
the size of the mask transition zone (bottom axis) and the mask 
roughness (top axis). 



ing from 0.93 D. A grey area in the plot underlines the range of 
Lyot stop sizes for which weak fluctuations of the coronagraph 
performance are observed. The CAPM coronagraph optimized 
with 90% LS is less sensitive to Lyot Stop reduction than that 
optimized with 100% LS. 

So, pupil undersizing weakly alters the CAPM coronagraph per- 
formance if this coronagraphic system has been optimized con- 
sidering the Lyot Stop size. Let us notice that in the case of a 
Lyot stop size larger than the geometric pupil diameter, the per- 
formance of the CAPM coronagraph decreases quickly as ex- 
pected since the rejected starlight is found there. 



6.4.2. Lyot stop misalignment 

Lyot stop misalignment has been identified as one of the 
issues for the coronagraphs on VLT-SPHERE. A misalign- 
ment less than 0.2% of the pupil diameter has been imposed 
to the selected cor onagraphs for the instrument requirements 
(Beuzit et al. 2006). In the following, we analyze the impact of 
the Lyot Stop misalignment with respect to the geometric pupil 
on the CAPM coronagraph performance. The LS size is taken 
here equal to 90% the relayed pupil diameter D. Once again, we 
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Fig. 23. Theoretical averaged intensity at 5 Ao/D achieved with 
the CAPM coronagraph as a function of the Lyot stop (LS) size. 
The dashed-line represents the edge of the geometrical relayed 
pupil. The grey area represents the tolerance range, defined in 
the text, for the CAPM coronagraph optimized with 90% LS. 



consider two CAPM coronagraphs: optimized with a 90% and 
100% LS. 

In Figure [24] the averaged intensity at 5 Aq/D provided by the 
DZPM coronagraph is represented as a function of the Lyot Stop 
misalignment. In the absence of pupil misalignment, the inten- 
sity level achieved by the CAPM coronagraph optimized with 
a 90% LS and 100% LS are about 4 TO" 8 and 8 TO -8 respec- 
tively. The result achieved by the CAPM coronagraph optimized 
to 90% LS is better than that optimized to 100% LS since we are 
facing a 90% LS in the coronagraphic scheme as mentionned 
above. 

As the pupil misalignment increases up to 0.05 D, the intensity 
level provided by the CAPM coronagraph does not really change 
significantly. Clearly, the CAPM coronagraph is not sensitive to 
pupil misalignment for a LS remaining completely inside the ge- 
ometric pupil. 

Beyond a 0.05 D pupil misalignment, the Lyot Stop starts to go 
outside the geometric pupil image where most of the rejected 
starlight is found. It results in a very high increase of the aver- 
aged intensity and an expected dramatic diminution in perfor- 
mance for the CAPM coronagraph. 

6.5. Star angular size 

A Sun-like star at lOpc is 1 milliarcsec across, corresponding to 
0.024 Aq/D and 0.125 Aq/D in H-band for a 8m-class telescope 
and the E-ELT respectively. Since phase mask coronagraphs are 
quite sensitive to star size dGuvon et al.ll2006h . we analyze the 
impact of this latter on the performance of our concept. 
In Figure 25, we plot the azimuthally averaged intensity profiles 
with the CAPM coronagraph for different star angular radii 9. 
All the profiles follow a a~ 3 law like the Airy diffraction pattern 
and the intensity values increase with 9, as expected. In Figure 
26 top and bottom plots, we represent the averaged intensity at 
5 Aq/D and IDA respectively as a function of the star angular 
radius. In the top plot, the intensity level goes beyond 10~ 7 for 
9 larger than 0.018 Aq/D. At larger angular radii, the intensity 
level increases as 9 2 . In bottom plot, a 3 Aq/D IDA is preserved 
for 9 lower than 0.025 Aq/D. At larger values, IDA evolves as 9. 
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Fig. 24. Theoretical averaged intensity at 5 Ao/D achieved with 
the CAPM coronagraph and a 90% Lyot stop, as a function of 
the Lyot stop misalignment. The dashed-line represents the edge 
of the geometrical relayed pupil. These results are obtained in 
the presence of a clear circular aperture. 

As a consequence, the finite size of the sun-like star at lOpc 
has little impact on the CAPM coronagraph when observed 
through an 8m-class telescope. However, for a 40m-class tele- 
scope, the finite size of this star clearly limits the coronagraphic 
performance . Thes e results confirm the limitations observed by 
Guyon et al. (2006) about the circular phase mask coronagraphs 
for the direct imaging of extrasolar terrestrial planets. In order to 
use the CAPM coronagraph at its optimal performance with an 
ELT, it is therefore more appropriate to observe more compact 
or more distant stellar objects. 

STELLAR RADIUS, 25% bandwidth 



0.00C» /D 




0.1 1.0 10.0 100.0 

Angular separation in X /D unit 

Fig. 25. Radial intensity profiles of the coronagraphic images 
theoretically achieved with the CAPM coronagraph for different 
stellar radii. 

7. Conclusion 

The CAPM coronagraph is an original concept composed of 
a colored pupil apodization, a DZPM in the following focal 
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Fig. 26. Theoretical averaged intensity at 5 Aq/D (top plot) and 
IDA (bottom plot) achieved with the CAPM coronagraph for dif- 
ferent stellar radii. 

plane and a Lyot Stop in the re-imaged pupil plane. This diffrac- 
tion suppression syst em constitutes an impr oved version of the 
DZPM coronagraph ( Soum mer et al.ll2~003bl) in which the grey 
apodization has been replaced by a colored one. This change 
allows us to increase the achromaticity of circular phase mask 
coronagraphs; indeed, a 2.5 mag contrast gain for the CAPM 
coronagraph is achieved with respect to the DZPM coronagraph. 
In addition, intensity levels of 2.2-10~ 7 and 2.210 s at 3 Aq/D 
and 5Aq/D respectively from the observed bright star are ex- 
pected with the CAPM coronagraph, in the case of a clear aper- 
ture and a 25% bandwidth. This performance situates the con- 
cept in an excellent position with respect to similar concepts re- 
viewed in the introduction to this paper. 

Manufacturing aspects of the colored apodizer have also been 
addressed, leading to a promising solution for the realization of 
our CAPM coronagraph using a nickel thin film for the colored 
amplitude apodizer. Optimizing the metallic layer thickness and 
mask parameters, we obtained the ultimate performance of the 
CAPM coronagraph. We consider to use a lens doublet combin- 
ing two glasses with the same refractive index, the same Abbe 
number and different anomalous dispersion to produce colored 
phase apodization. An encouraging solution has been found for 
the visible range and studies are on-going to address other spec- 
tral bands with the same concept. 
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A sensitivity analysis of the CAPM coronagraph to different er- 
rors has also been performed, leading to the determination of the 
error level which can be accepted before the performance of the 
device is altered. We find in particular that for a given amount 
of rms wavefront errors, the concept is less sensitive to low- 
order than to high-order aberrations. While an intensity level at 
5 Ao/D of 1(T 7 is ensured for about 6 • 10~ 3 /to of tilt or defocus, 
the same intensity level requires spherical aberration less than 
3 • KT 3 i ( ). In the case of random aberrations with a v 2 PSD 
distribution, this intensity level requires wavefront errors below 
10~ 3 Ao. Similar observations are made in the case of close-in 
imaging capabilities, quantified by the Inner Detectability Angle 
(IDA). 

Constrains on DZPM design are quite relaxed at level of mask 
roughness (76 nm rms in the visible spectral range). However, 
they can prove to be critical in the case of the mask transi- 
tion width where a value better than 1 fim is required if one 
wants to use our CAPM coronagraph within an exoplanet im- 
ager with a reasonable focal ratio (F less than 100) in the visible. 
Technological improvements of mask manufacturing are under 
development and some first prototypes are expected to be tested 
in our laboratory soon. In addition, we show that an optimization 
of the CAPM coronagraph with a 90% Lyot stop is helpful for 
this device to be less sensitive to pupil stop misalignments. 
Furthermore, the averaged intensity at 5 Ao/D remains lower 
than 10~ 7 with our concept for a stellar angular size smaller than 
0.018 Ao/D. 

Central obstruction and spiders effects have not been addressed 
here. They will be presented in a future paper with a design 
study and performance analysis of the CAPM coronagraph in 
the framework of the E-ELT. 

Finally, the CAPM coronagraph constitutes a very attractive con- 
cept for the planet imagers of the next decade thanks to its 
small inner working angle, its great ability to suppress broad- 
band starlight and the relative simplicity of its design. The next 
step will now consist in manufacturing all the components of the 
CAPM coronagraph and performing an experimental demonstra- 
tion of the concept. We are currently preparing a visible corona- 
graph testbed at LAM in Marseilles for this purpose. 
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